Structure of the most singular vortices in fully developed turbulence 
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Using high Reynolds number experimental data, we search for most dissipative, most 
intense vortices. These structures possess a scaling predicted by log-Poisson model for the 
dissipation field e r . These new experimental data suggest that the most intense structures 
have co-dimension less than 2. The log-Poisson statistics is compared with log-binomial 
which follows from the random /3-model. 
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It is known, at least from numerical simulations, that 
the large-amplitude dissipation occurs around vortex 
tubes in turbulence. We thus expect some structure to 
exist in a signal that characterizes the large values of 
the dissipation field. There are some statistics, although 
very incomplete ones, on the distance between the vor- 
tex tubes, the size of vortex tubes, etc.. The largest 
value of the dissipation is also important in determin- 
ing the resolution of DNS [1] . It would be therefore of 
interest to provide a direct experimental study of the dis- 
sipative field extremal values. On the other hand, very 
large values of the dissipation field correspond to inter- 
mittency. Traditionally, the latter is expressed through 
so-called intermittency corrections to the exponents for 
the structure functions, (\u(x + r) — u(x)\ p ) ~ r^ 35 , where 
u is the longitudinal velocity, and ( p = p/3, [2]. Thus, 
these corrections result in ( p — p/3 + r(p/3). A theory, 
incorporating the intermittency, the refined similarity hy- 
pothesis [3], links the statistic of these corrections with 
the statistic of the dissipation field e r , the energy dissipa- 
tion averaged over a ball of size r. Namely, (e£) ~ r r ( p \ 
Many models have been proposed to explain intermit- 
tency. It was originally suggested that the statistics of 
e r is log-normal [3]. More recently, She and Leveque [4] 
(hereafter SL, see also [5], [6], and recent study [7]) have 
proposed log-Poisson statistics for the dissipation field, 
with agreement with the experimentally found ( p in [8], 
[9]. These experimental exponents are obtained in Ex- 
tended Self-Similarity approach, which is useful because 
of extended scaling range. 

The simplest idea to study large values of dissipation 
is to measure the maxima. However, for many distri- 
butions a maximum of a big array can be "anything" , 
or arbitrary large. This is true, for example, for Gaus- 
sian statistics. The same is true for log-normal distribu- 
tion. To see this, recall that, first, t(j>) = —d p (p — 1), 
d p = D — D p , where D p are so-called generalized di- 
mensions [10]. Second, studying maxima is in a way 
equivalent to measuring asymptotically high moments, 
doo = linip^oo {— t(p)/(j) — 1)}. For the log-normal dis- 
tribution, r(p) = —/j,/2p(p — 1), and therefore asymp- 



totically, d p = (fi/2)p — > oo. Remarkably, the Poisson 
statistics provide some distinctive maximum. To see this, 
recall that for the Poisson distribution (see, e.g., [11]), 



>£7a!, a = 0,1, 2,. 



aa+b 
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and a < 0, the maximal value is defined through b = 
lnmaxe r . In order to specify b, we calculate the mo- 
ments, (e^). Noting that (e r ) = 1, we get b — £(1 — e a ), 
where £ = Cln(£/r), I being external scale, and C is a 
constant. As a result (of calculation of the moments), we 
get, 



r(p) = c[i - (i - j/cy 



PI, 



(2) 



where 7 = C(l — e a ) [12]. Using (2), it is easy to show 
that this time doo = 7, which is a finite number. 

SL is recovered from (2) if C = 2, 7 = 2/3, so that a = 
In (2/3), C being the co-dimension of most dissipative 
structures, and 7 is defined by the dissipation rate, i.e., 



inverse time-scale, 1 jt r 



-2/3 



[4]. 



The meaning of C becomes even more clear directly 
from (1): the most intense fluctuations correspond to 
a = (as a < 0), so that the probability P(a = 0) = 
= (rji) c = (r /£y D ~ H °\ D is dimension of space 
(— 3). Thus the Hausdorff dimension for most dissipative 
structures in SL theory, Hq = 1, i.e., the structures are 
filaments. On the other hand, using expressions for b, £ 
and 7, we now rewrite (1) as follows, 



(3) 



Putting a = in (3), we can see that the most intense 
structures are expected to scale ~ r~ 7 . Thus, for the 
the log-Poisson statistics, the maxima of e r (x) are not 
"anything" , and they are supposed to be self-similar. 

This scaling is proved to be possible to verify exper- 
imentally. We used 10 million points of atmospheric 
data, with an estimated Taylor microscalc Reynolds num- 
ber 9540, (experiment A) and 40 million points for both 
longitudinal and transfer velocities (experiment B). The 
data are treated in spirit of Taylor hypothesis, that is, 
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the time series is treated as one-dimensional cut of the 
process. The dissipation rate can be written as 



e(x) = u(diVjdiVj +diVjdjVi), 



(4) 



(summation over repeating induces). The second term 
on the rhs vanishes after averaging, for homogeneous in- 
compressible turbulence. The first term consists of 3 lon- 
gitudinal and 6 transverse components. Therefore, it is 
natural to present the dissipation as 



e c (x) = 3 (d x v x (x)f + 6 (d x v y (x)) 2 



(5) 



For isotropic turbulence, {(d x v y (x)) 2 ) = 2((d x v x (x)) 2 ) , 
and therefore, following [13], we may consider three types 
of dissipation, longitudinal, transverse, 

£i(x) = 15z/ {d x v x {x)f , e t {x) = (15/2)i/ {d x Vy(x)f , 

(6) 



and combined, (4). 



than 2/3. In order to compare with a "regular" ran- 
dom process we generated a Gaussian process ui g with 
correlation function coinciding with experimental, i.e., 
(bj g {x + r)uj g (x)) = (d x v x (x + r)d x v x (x)). Then, the "dis- 



sipation" e (f) = u 2 g , and e [ r 9} = l/r f^/i U (a) {x')dx'. 
Corresponding calculation for the maxima are reported 
in Fig. 1. If any scaling can be extracted from the Gaus- 
sian process, it would be at large asymptotic distances, 
and the scaling is trivial, 7 = 0, meaning no singularity. 

Figure 2 presents the scaling for experiment B. This 
time, the scaling holds for almost 6 decades. The scaling 
exponent for longitudinal dissipation is again 0.61 ±0.01, 
while for the combined dissipation (which is quite close to 
the transversal dissipation) the exponent is 0.57 ± 0.01. 
In both Figs. 1 and 2 there is also A-scale. Note that 
there is characteristic transfer region at r/\ « 1. We 
may interpret it as a transition to the inertial range [14] , 
which is formed due to the fact that the vortices are 
expected to have scales between Kolmogorov microscale 
rj and Taylor microscale A. 
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FIG. 1. Scaling for most intense structures. The power law 
fitting of the experimental data (solid thick line) has been 
extended to reach unity (solid line), where it is supposed to 
match with SL scaling. The distances are given in terms of 
Kolmogorov micro-scale rj, and in units of A. 

We will deal with coarse-grain dimensionless dissipa- 
tion, 



1 [ x+r/2 e(x')dx' 
r J x -r/2 (e) 



(7) 



and maxima of e r can be measured. Note that this mea- 
surement is meaningful because (7) contains some aver- 
age. Figure 1 shows longitudinal scaling for the experi- 
ment A, which holds for 4.5 decades. The deviation from 
SL is small, and we recall that SL suggest that there is 
no anomalous scaling for t r . This small deviation in Fig. 
1 can be interpreted as anomalous persistence of the ed- 
dies, which is indeed observed [15], see also discussion in 
[16]. The value of 7 is 0.61 ± 0.01, only slightly smaller 
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FIG. 2. Experiment B scaling for longitudinal and 
transversal dissipation, including combined dissipation. 

Although the experimental 7 is not that different from 
2/3, the value of the other parameter C is quite sensi- 
tive to that difference. In order to find C we substitute 
7 from our measurements into (2), and use computer 
routines to find a best fit for these data with free pa- 
rameter C and the exponents C,p ESS ^ from experiment 
[8], [9]. We start with the longitudinal dissipation (6) 
from experiments A and B (recall that 7 is the same for 
them). As a result, we find C = 1.67 and a = —0.45 (cf. 
In {2/3} = —0.41). With these parameters, the deviation 
of these calculated exponents ( p 

tal exponents 
pare: for SL, 



(e) 



from the experimen- 
ted ~(p ESS) ) 2 ) = 0.0063. To com- 



UC 



(SL) 

p 



AESS) 



) 2 ) = 0.0078. The Q 

(SL) 



(e) 



exponents seem to be "better" than Q , but consider- 
ing that the experimental exponents have errors of about 
±1% [9], we conclude that these exponents are similar. 
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Note that if we substitute in (2) the value of 7 = 0.61 



and put C 
we have, \l {{Q 



2, then, for the obtained exponents, Q, 



(C=2) 



(C=2) 



( { p ESS) )2) = 0.050, much too high. 
Consider now the combined dissipation, defined in (5). 
The best fitting with this parameter fixed results in 
C = 1.43, and a = —0.50. This time, the deviation of 
the computer generated spectrum from the experimental 
is = 0.0071. Figure 3 shows r(p/3) from experiment, and 
for different theories. It can be seen that all the curves 
collapse into one, corresponding to the experiment, ex- 
cept that one with 7 from our measurements, and C = 2. 
This illustrates that the data are indeed sensitive to the 
measured 7, that is to its (small) difference from 2/3. 
The codimension C = 1.43 corresponds to Hq = 1.57. 
This value of H a > 1 seems to be consistent with the 
distinction between persistent vortical filaments and the 
dissipative structures associated with regions of strong 
strain [18]. That means that the most dissipative struc- 
tures consist not only of filaments, but in part of sheets, 
or filaments convoluted into complex structures, cover- 
ing more than 1 dimension. According to the intersec- 
tion theorem [19], that D - H = £> (m) - H^ m \ where 
£j(m) j g di mens i on f the measurements (in our case 
£)(m) _ an( j _ corresponding measured Haus- 

dorff dimension. It is clear from this formula that, if 
H < 2 then i^ m) < 0. This actually means that the 
dimension H < 2 cannot be detected in ID measure- 
ments directly, and therefore our conclusion is inevitably 
indirect. Indeed, it is obtained from spectrum (2), really 
formed in 3D, but projected into ID assuming isotropy. 
Therefore, it would be important to measure the Haus- 
dorff dimension in 3D simulations directly Another rea- 
son for that is the surrogacy issue [20] . 

These statements about the dimensions of the most in- 
tense structures can be also formulated for log-binomial 
distribution, and, it is known that the Poisson process 
is a limit of the binomial distribution for "rare events" . 
In particular, the Poisson distribution can be obtained 
from the random /3-model [17] by a suitable limiting 
process [5], [21]. Let (3 take two values, W — (3\ with 
probability x, and W = [3 2 with probability 1 — x, and 
(3\x + /?2(1 — x) = 1 (in order to have t(1) = 0). Let also 
Pi < 1 < 02- Then, on the n-th level, the distribution is 
binomial, that is, W n = e n = (3™ (3% ~ m with probability 
(™ )x m {\ - x) n - m . Hence, (e£) = [x/3{ + (1 - x)0%] n . 
Taking into account that n = In (r/^)/lnr, T being the 
ratio of successive scales, we obtain, 



t( P ) = 



lnfcffi + (l-:c)/%] 
lnT 



(8) 



In [5] and [21], T was treated as a free parameter. It was 
shown that, if T = 1 — x/C, (3\ = 1 — 7/C and x — > 0, 
then (3 2 ~ 1 + x-f/C, and (8) reduces to (2). 
The most intense structures on n-th level, 



/3 2 n 



G) 



ln/3 2 /lnr 



-7/3 



(9) 



cf. (3). On the other hand, the probability of these 



maxima, 



p = {i-xY = ^ 



r \ ln(l-x)/lnr 



-0 



(10) 



In particular, if T = 1 — x/2 and x — > 0, then Cp = 2 
[22]. 

If we do not treat T as a free parameter, and consider 
that it is a fixed number, then the log-binomial distri- 
bution generally cannot be reduced to the log-Poisson 
PDF: in particular, if x is small, then, r(p) ~ x — > 0, 
and thus the intermittency is negligible . As in our case 
the division level n ^> 1, the log-binomial distribution be- 
comes essentially log-normal with maximum at m = xn 
[11], [25]. However, the log- normal distribution has many 
shortcomings, and it has been repeatedly criticized when 
used to explain anomalous spectrum [26], [27]. Neverthe- 
less, the spectrum (8) does not even look like log-normal 
(for which r(p) = —(fi/2)p(p — 1)) and rather behaves 
like log-Poisson for p 3> 1. Indeed, according to (8), for 
P» 1, 



T{p)=C p + C^-vrip, C\ 



(1-x) lnT 



<0, (11) 



(3 — (3\/ @2- This spectrum resembles (2); and the con- 
stants in (11) happen to be numerically close to corre- 
sponding numbers in (2). The reason for such a dramatic 
difference with log-normal distribution is as follows. For 
binomial distribution, 



m=0 



(12) 



S — (3 P <C 1 for large p. Then S m decreases dramati- 
cally with increasing to, and therefore the terms of the 
sum (12) of maximal probability, at to ~ xn, where nor- 
mal distribution if formed, do not contribute substan- 
tially. In contrast, only the first few terms of this sum 
(responsible for "rare" and very intense events) really 
contribute. Thus, effectively, the distribution works like 
a Poisson distribution. To see this explicitly, consider 
a probability distribution (™ )x m (l - x) n ~ m 5™/A, where 
A is a normalization constant, A = (xSq + 1 — x) n , and 
^0 = (/3i//?2) Pt \ Po 3* 1- Then, for large n we express 
the factorials entering the binomial coefficients through 
Stirling formula (except for to!, because to is not neces- 
sarily large), to get, 



Pi(m) = jUx m (l-x] 



n—m cm 



-£o50 (13) 



where £0 = n5ox/{l — x). For p > po, the sum (12) can 
be written as A(3 2 p £m=o Pi(m)S m( - p - po \ and thus the 
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distribution effectively corresponds to the Poisson distri- 
bution. 

Let us take the random j3 model "for real", that is, 
consider the Poisson distribution as an approximation to 
the binomial, as in (13). Then, we may consider that the 
model is realized as follows. Denote the number of divi- 
sions of each volume by N. Then T = l/N 1 ' . We now 
multiply the values of m divisions by (3\ and multiply the 
remaining N — m divisions by /?2- This actually means 
that the probability x = m/N, and 1 — x = (N — m)/N). 
A particular case of N = 2 and D = 1, i.e., x = 1/2, cor- 
responds to the model proposed in [23] . Then, according 
to (10), P = (r/£) D , i.e., the Hausdorff dimension Hq 
is = 0, while 73, defined from (9), = L>ln/3 2 /m2. The 
case Pi = returns us to the /3-model [24]. In that case, 
j3 2 = 1/(1 - x), and 7j8 = Cp = Din (1 - x)/ In (1/2). 



intense structures geometry in fully developed turbu- 
lence. The PDF's of the exponents of the dissipation 
field are compared with the log-Poisson distribution to 
show a good agreement with the theory. The log-Poisson 
statistic can be considered as a limiting case for the log- 
binomial distribution appearing in random /3-model. We 
estimated the parameters of the log-binomial distribution 
with 7 found in our measurements, and to fit the expo- 
nents for the structure functions found elsewhere. We 
conclude that the estimated Hausdorff co-dimension of 
the most intense structures is less than 2. 

We appreciate numerous comments made by S. 
Boldyrev, Z. Mikic, and R. Rosner. 
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FIG. 3. Intermittency corrections from experiment [9], and 
from other theoretical models. 

In general case, we may take D = 3, and we are deal- 
ing with (8) with 7^ given from our measurements (so 
that /?2 is defined according to (9)). We thus are left 
with two free parameters, N, an integer, and x = m/N, 
where m is also an integer. These numbers can be found 
with help of computer search to fit experimental data 
[8], [9] in optimal way. As a result of this search, we 
get: For 7 = 0.57 (combined longitudinal and transver- 
sal dissipation), N = 11, x = 7/11. With these pa- 
rameters, the deviation of the spectrum from experimen- 
tal is 0.0098, quite satisfactory. Indeed, the correspond- 
ing i~(p/3) depicted in Fig. 3 is indistinguishable from 
other approximations which collapse to the experimental 
data. For 7 = 0.61 (longitudinal dissipation), N = 4, 
x = 2/4 = 1/2, Cp = 1.5, the deviation is 0.0115, still 
okay. As mentioned, at p ^> 1, the log-binomial spec- 
trum (8) is essentially reduced to the log-Poisson spec- 
trum (2), and therefore we would prefer to consider the 
log-binomial distribution to be more general. 

In conclusion, one of the predictions of SL theory about 
the scaling of maxima ~ r~ 7 is experimentally confirmed. 
This makes it possible to make a better estimate of the 
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